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Abstract 

We investigate the quantum versions of a one-dimensional random 
walk, whose corresponding Markov Chain is of order 2. This corresponds 
to the walk having a memory of up to two previous steps. We derive the 
amplitudes and probabilities for these walks, and point out how they differ 
from both classical random walks, and quantum walks without memory. 

1 Introduction 

"Standard" One Dimensional Discrete Quantum Walks (also known as Quantum 
Markov Chains) take place on the state space spanned by vectors 

\n,p) (1.1) 

where n E Z (the integers) and p € {0, 1} is a boolean variable (see [HHITO] for a 
comprehensive treatment). The second variable p is often called the 'coin' state 
or the chirality, with representing spin up and 1 representing spin down. It 
is the quantum part of the walk, while n is the classical part. One step of the 
walk is given by the transitions 

\n,0) — > a\n- 1,0) +&|ra + l,l) (1.2) 

\n,l) — > c\n - 1,0) + d\n + 1,1) (1.3) 

where 

(c !0 e5f/(2) ' (L4) 

the group of 2 x 2 unitary matrices of determinant 1. These walks have been 
well studied, and their asymptotic behaviour well analyzed [Ml [HI EJ [21 [1] . 

The corresponding classical walk is represented by a Markov Chain whose 
transition matrix is tridiagonal, with zeroes along the diagonal, and 1/2 along 
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off-diagonal (for the fair coin): 
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(1.5) 



In this paper we investigate quantum walks "with memory" : The state space 
is spanned by vectors of the form 



\n r ,n r -i,...,n 2 ,rii,p) 



(1.6) 



where rij = nj-i ± 1, since the walk only takes one step right or left at each 
time interval, nj is the position of the walk at time t—j + 1 (so n\ is the current 
position). The transitions are of the form 

|n r ,n r _i, . . . ,n 2 ,n 1 ,0) — > a |n r _i, ... ,n 2 ,n 1 ,n 1 ± 1,0) 

+ 6|n r _i,...,n 2 ,ni,ni ± 1, 1) (1.7) 

|n r ,n r _i, . ..,n 2 ,ni, 1) — > c|n r _i, . . . , n 2 , m, rii ± 1, 0) 

+ d\n r -i,...,n 2 ,ni,ni ± 1, 1) (1.8) 

In analogy with the definition for Markov Chains, we call r the order of the 
quantum walk. 



2 Order 2 walks 

The state space is composed of the families of vectors 

|n— 0) , \n — l,n, 1) , |n + l,n, 0), |n + l,n, 1) 



(2.1) 



for n G Z. In what follows, we will refer, for obvious reasons, to \n — 1, n,p) as 
a right-mover, and to \n + l,n,p) as a left-mover. Following pQ, it will suit us 
also to split the transitions (Eq. 1.7 1.8 1 into two steps, a "coin flip" operator 
C and a "shift" operator S: 



C 

c 
s 



\n 2 ,ni,0) — > a\n 2 ,ni,0) +b\n 2l n 1 ,l) 
|n2,rci,l) — > c|n 2 ,rii,0) +d\n 2 ,n 1 ,l) 
l«2,«i,p) — > \nx,ni ± 



(2.2) 
(2.3) 
(2.4) 



We investigate in what follows the possibilities for the shift operator 5*. Suppose 
S sends both \n — 1, n, 0) and \n + 1, n, 0) to the same vector, say |n, n + 1, 0) 
(thus, in our parlance, for p = 0, it sends both left and right movers to right 
movers). One observes immediately that this is not really a 2nd. order chain 
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(its behaviour does not depend on 712, only on rii). Indeed, on our state space 
\ri2, Ti\ , p) it is not even unitary (even though it would be on the state space \n, p) 
of an order 1 walk). For the behaviour with p = 1, we have two possibilities: 

1. S sends both \n — 1, n, 1) and \n + 1, n, 1) in the same direction (whether 
left or right). In this case, again S behaves as a first order transition, and 
the whole analysis is that of a 1st. order quantum walk. 

2. S sends \n — 1, n, 1) and \n + 1, n, 1) in different directions. So, for p = 1, 
S behaves like a 2nd. order chain. In this case, it turns out that the com- 
bined behaviour does not give an invertible transition: i.e. the transition 
matrix is not unitary. 

Because of these arguments, to construct a bona fide 2nd. order walk, S 
needs to send \n — 1, n,p) to a different state than it sends \n + 1, n,p), for both 
values of p. The four possibilities are described in Table [TJ There is a simple 



Initial State 


Final State 




Case a 


Case b 


Case c 


Case d 


|n- l,n,0) 


\n,n + 1,0) 


|n,n + 1,0) 


\n, n — 1, 0) 


|n, n — 1, 0) 


\n- l,n,l) 


n, n + 1, 1) 


|n,n- 1,1) 


n, n + 1, 1) 


|n,n- 1,1) 


|n + l,n,0) 


|n, n — 1, 0) 


|n, n — 1, 0) 


|n,n+ 1,0) 


|n,n + l,0) 


\n + l,n, 1) 


\n,n- 1,1) 


|n, n + 1, 1) 


\n,n- 1, 1) 


|n, n + 1, 1) 



Table 1: Action of shift operator 5" 



way to view these cases, as follows. Depending on the value of the coin state p, 
one either transmits or reflects the walk: 

Transmission corresponds to \n — 1, n, p) — > \n,n + l,p) and | n + l,n,p) — s- 
\n,n— l,p) (i.e. the particle keeps walking in the same direction it was 
going in) 

Reflection corresponds to \n—l,n,p) — > \n, n— l,p) and \n + 1, n, p) — > 
\n, n + l,p) (i.e. the particle changes direction) 

We re-phrase in Table [2] the action of 5* described in Table [l] 



Value of p 


Action of S 




Case a 


Case b 


Case c 


Case d 





Transmit 


Transmit 


Reflect 


Reflect 


1 


Transmit 


Reflect 


Transmit 


Reflect 



Table 2: Action of shift operator S 



2.1 Initial Conditions 

We must clarify how to initialize the walk, since at the very beginning, we 
cannot run a 2nd. order chain without any history. "Starting" at position - 
1, we then move to position (which can be done using a first order quantum 
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walk). This creates the state | — 1, 0, 0), and from there on we can run the second 
order operations described above. 



2.2 The Hadamard Walk 

We observe that Cases (a) and (d) do not lead to any interesting features. In 
Case (a), the particle just moves uniformly right or left, depending on the initial 
state. If the initial state is a superposition of left- and right- movers, the walk 
progresses simultaneously right and left. For Case (d), the walk "stays put", 
oscillating forever between n and n + l for some value of n. In both these cases in 
fact, the coin flip operator C plays no role (since the action of S is independent 
of p), so there is nothing quantum about these walks. 

However, cases (b) and (c) do yield results of interest. To analyze these, we 
choose a particular coin flip operator C corresponding to the Hadamard walk: 

Classically C sends |n 2 , ni,p) to either |n 2 , n ij 0) or |n 2 , n%, 1) with equal prob- 
ability 1/2 (fair coin toss). 

Quantumly 

C: K, ni,0) — > -^=(|n 2 ,ni,0) + \n 2 ,n u l)) (2.5) 
C : |Ti2,ni,l) — ► -=(|n 2 ,wi,0) - \n 2 ,ni, 1)) (2.6) 



The equations 2.5 and 2.6 correspond to a = b = c = —d = l/y/2 which is 
known as the Hadamard walk. 

For both cases (b) and (c) it should be clear that in the classical case, we 
end up with the standard (classical) random walk: In each case, transmission 
and reflection just correspond to picking one of two different choices (right or 
left) at each step. 

Let us consider case (c): The first few steps of a standard quantum (Hadamard) 
walk starting at position n would be 

1 



,0) 



V2 



(|n-l,0) + |n + l,l» 



1 



2^2 



-2,0) + \n, 1) + |n, 0) - \n + 2, 1)) — > 

(|n - 3, 0) + \n - 1, 1) + |n - 1, 0) - |n + 1, 1) 



(2.7) 
(2.8) 



+ \n - 1, 0) + \n+ 1, 1) - \n+ 1, 0) + \n + 3, 1)). (2.9) 

Thus after the third step of the walk we see destructive interference (cancellation 
of 4th. and 6th. terms in expression 2.9) and constructive interference (addition 
of 3rd. and 5th. terms in expression 2.9). However for case (c) the first few steps 
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are for example 



|n-l,n,0) — -=(|n,n-l l 0) + |n,n + l,l)) — > (2.10) 
v 2 

-(|n- l,n,0) + | n -l,n-2,l> + |n + l,n,0) - |n + l,n + 2,l)) — > (2.11) 
|n,n- 1,0) + |n,ra + l,l) + |n-2,n- 1,0) - \n-2,n- 3, 1) 



2a/2 



+ |n, n + 1, 0) + |n, n - 1, 1) - |n + 2, n + 1, 0) + |n + 2, n + 3, 1)) — ► 

(2.12) 

^(|n- l,n,0) + |n- l,n- 2,1) + |n + l,n,0) - |n + l,n + 2, 1) 

+ |n - 1, n - 2, 0) + \n - 1, n, 1) - \n - 3, n - 2, 0) + |n - 3,n - 4, 1) 

+ |n + 1, n, 0) + |n + 1, n + 2, 1) + \n - 1, n, 0) - |n - l,n - 2, 1) 

- |n + 1, n + 2, 0) - \n + 1, n, 1) + \n + 3, n + 2, 0) - |n + 3, n + 4, 1)). 

(2.13) 

After three steps, there is no interference (constructive or destructive), but the 
interference appears after step four (e.g., in expression 2.13 we can cancel the 
2nd. and 12th. terms, and we can add term 3 and term 9, etc.). Thus we can 
see this walk differs both from the classical random walk and from the standard 
(Hadamard) quantum walk. 



3 Amplitudes 

(This section follows closely the approach taken in Appendix A of [3]). We now 
derive analytical expressions for the wavefunction amplitudes in case (c) of table 
[2j using as quantum coin flip the Hadamard transition |2.5| and |2.6| 

For the 1-dimensional walk, we view the progression as a sequence of left 
(L) and right (i?) moves. In general there are many paths to reach a particular 
final state: We need to sum over the amplitudes of these different paths (with 
appropriate phases) to obtain the amplitude for that final state. 

As a quick example, for the classical case, what is the probability of ending 
at position 1 in a 3-step walk that starts at the origin? The possible walks 
ending at 1 are LRR, RLR or RRL. The total number of possible 3-step walks 
is 2 3 = 8. So the probability of finishing at positions 1 is 3/8. 



In the notation of expression 2.1 let our initial state be | — 1,0,0) (so the 
walk starts at the origin) and let us take n steps in the walk. It should first of 
all be obvoius that as in the classical case, if n is odd/even, we can only finish 
up at an odd/even integer position (respectively) on the 1-dim lattice. Let Nl 
be the number of left moves, and Nr the number of right moves. 

Lemma 3.1 We refer to an 'isolated' L (respectively R) as one which is not 
bordered on either side by another L (respectively R). Let (respectively N^) 
be the number of isolated Ls (respectively isolated Rs) in the sequence of steps 
of the walk. Then, the quantum phase associated with this sequence is 
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Proof In what follows, we first analyze the sequence of Ls (identical arguments 
will apply to the Rs). An isolated L does not contribute to the phase, nor 
does the pair LL bordered by Rs. The first sequence of Ls that can contribute 
is LLL: In our previous language, this corresponds to transmit followed by 
transmit. After the first L, the coin state is 0, after the second it is 1, and after 
the third it is 1. It is the transition from 1 to 1 in the coin state that gives the 
factor of —1 from the Hadamard walk. 

So in general, a sequence of j LL . . . Ls will give a phase contribution of 
(-1)3 for j > 2. 

Now examine clusters of Ls of size greater than 2. If we have 2 such clusters, 
we can move one L from the first cluster to the second, without changing the 
overall phase contribution. In such a move, the contribution of the 1st. cluster 
decreases by a factor of -1, while that of the 2nd. increases by the same factor. 
Suppose we repeat this process, to shrink all but one of the large clusters to 
clusters of size 2. We end up with a sequence that looks like 



. . . RLR . . . RLR . . . RLLR . . . RLR . . . RLLR ...R LLLLL ...L R... 

S v ' 

One large cluster of Ls 

(3.2) 

Denote by Cl the total number of L clusters. Then the total number of L 
clusters of size 2 is Cl — N\ — 1. So, the size of the one large cluster of Ls is 
N L -Nl- 2{C L - N X L - 1) = N L + Nl - 2C L + 2. Its phase contribution is 
therefore (-1)^+^ . 

Since analogous arguments apply for sequences of Rs, the total phase con- 
tribution is (-1)Nl+n r +nI+n 1 r _ i 

After an n- step walk, we want to know what is the probability the particle 
is in position k. From previous arguments, (—1)" = (— l) k and — n < k < 
n. Four possible final quantum states correspond in our model to the particle 
terminating at k: 



\k-l,k,0) , \k-l,k,l) , \k + l,k,0) , |fc + l,fc,l) (3.3) 

sequence sequence sequence sequence 

ending . . . LR ending . . . RR ending . . . RL ending . . . LL 

Let us denote by a^^R, a^RR, OLkRL, o-hll the amplitudes of these 4 states in the 
final wavcfunction ^ . Then the probability when we measure of finding the 
particle at position k is 

Wlr\ 2 + Kfli?| 2 + Wrl\ 2 + Kll| 2 (3.4) 

Before calculating the amplitudes, we need another technical lemma. 

Lemma 3.2 Consider a composition (ordered partition) of the integer n into 
C parts, and let N 1 be the number of ls in the composition. Then either 

1. n = C = N 1 
or 
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2. max(0, 2C - n) < N 1 < C - 1. 

Proof Case 1. is trivial: It is the composition of n into N 1 Is. For case 2., the 
upper limit is also trivial: The largest number of individual Is we can get is 
C — 1, which is the composition 

n= 1 + 1 + 1 + ■•• + !+(» -((7-1)) (3.5) 

s V ^ 

(C— l)tcrms 

For the lower limit, assume C < n/2. It is always possible to write down a 
composition with few terms, without using any Is. Specifically, we can write 
the first C — 1 terms as 2, and the last term as the remainder (n — 2(C — 1)), 
which is greater than 2 by assumption. 

Now assume C > n/2. The least number of Is in the composition is obtained 
by writing as many 2s as possible. Suppose we have r 2s, and the other terms 
are 1. Then 2r + iV 1 = n. Since r = C- N 1 , we have that N 1 = 2C - n, and 
the result follows. | 

We define the combinatorial symbol 



a 



, ,,awc-a-l^ a! (c-a - 1)! 

" * 1 6 j M(a-6)!(c-2a + &)!(a-6-l)! [ ' 

Theorem 3.3 The amplitudes akLL, o-kLR, a kRLi&kRR for the final states given 
in Equation \3.3\ are 



Nr 

+N} N j(N R -Nj + 1) 
N R (N R + 1) 

Ni = „,x(l, 
2KT R -S L +2) 





(3.7) 



N L -l O-l C-l 



2 f ^=E E E i-ir +N t +N » 



C = 2 si = max(l, N R = max(l, 

2C-N,) 2C-N„) 




c 

K. ) +s Nl , Nr 



i AT 1 / Nr \ Nl - 1 lJV i / N L 

E ^ + E (-l) NR+N «^[ *r 



21V R -JV L ) 2N L -N R ) 

(3.8) 
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N L -1 C-l C-2 




Nl(Nl - 1) 



2N H -JV I/ +2) 

N L -2 



N R = 

2JV i -N R -2) 





+ E (-1) W M (3.9) 



ATr — 1 C-l C-l 



2 f ™=E E E (-ir +N t +N « 

C=l JV^ = max(l, N]j = max(0, 



Nl{C-Ni) 




Nl = max(0, 




where k — Nr — Nl , ?i = Nr + Nl — 2 and 5 is the standard kronecker delta 
function (5 Ptq = 1 ifp = q and zero otherwise) . 

Proof Because of its slightly lengthy and technical nature, we relegate the proof 
to Appendix A. | 



4 Simulations and Analysis 

We show in Figures [l] and [2] the amplitudes for the 3 different kinds of walks 
(classical, quantum, quantum with memory). The simulations are carried out 
using AXIOM [7]. For completeness, we include in Appendix B the commented 
code for generating the Quantum Walk with memory. 

In figures [l] and [2j the initial states for the three cases are |0), |0,0) and 
| — 1,0, 0) (by abuse of notation, the ket vector here |0) represents the classical 
case). As has been pointed out by a number of authors (see e.g. [HI Appendix 
A]) in the quantum case we can choose a more symmetric initial state (still 
of course representing the particle starting at the origin). In general this will 
give rise to a different probability distribution. For the quantum walk we start 
at ( 1 , 0) + 1 0, 1) ) / a/2 and for our walk with memory, we start at (| — 1,0, 0) + 
| — 1,0, 1) + 1 1, 0, 0) + 1 1 , 0, l))/2. The probability distributions for these cases 
(for a 40-step walk) are plotted in figure [3] 

What is immediately noticeable is the high probability that the quantum 
walk with memory stays at the origin (even after 40 steps, it has more than 
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Figure 1: Probability Distribution 
after 10 steps 



Figure 2: Probability Distribution 
after 40 steps 
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Figure 3: Probability Distribution after 40 steps for symmetric initial state 



50% chance of being found at the origin) . In the terminology of Konno [TT] , we 
say the particle is localized at the origin. Also of note are the smaller peaks 
that occur quite a distance from the origin (at ±6 in Figure [T] and at ±28 in 
Figure [2]). The distribution is symmetric about zero, except for the one specific 
case Nft = ± 1 (i.e. in a walk with an odd number of steps, the probability 
of finding the particle at positions ±1 is not the same). 

Claim As the quantum walk with memory becomes infinitely long, for even n 
there is still a chance of over 50% of finding the particle at the origin! 
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Proof The proof proceeds by setting n = 2j and using an inductive argument 
on j (the particle can only be at the origin for an even number of steps) . Let us 
denote by a,k**(n) the dependence of the amplitude on the number of steps n, 
where ** is one of LL, LR, RL, RR. The argument focuses on the dependence of 
S = {doLR(n),a RL(n)} on their equivalents two steps earlier = {a aLR (n — 
2),a QRL (n-2)}. 

Base Case For n = 2, a LR,(2) = aofli,(2) = 1/2 are the only terms contribut- 
ing to the probability at the origin. 

Inductive Step Let us consider aoLR(n) and ao.RL(n). Assume the amplitudes 
a-0LR.{n — 2) and oorl{ti — 2) are both positive and sum to 1 (as in the 
base case). 

Amplitude of | — 1,0,0) This corresponds to a L R {n). There are 2 con- 
tributions from ao**(n — 2): 

Contribution from aoL R (n — 2) The particle moves left and then 

right. The phase contribution stays positive. The amplitude 

factor is (l/\/2) 2 = 0.5. 
Contribution from a 0R L(n — 2) The particle moves left and then 

right. The phase contribution stays positive. The amplitude 

factor is (l/\/2) 2 = 0.5. 

Thus the total amplitude contribution is 0.5a LR( n —2)+0.5a aRL (n — 
2) = 0.5a OL fl(n - 2) + 0.5(1 - a QLR (n - 2)) = 0.5 

Amplitude of |1,0, 0) This corresponds to a 0R L(n). There are 2 contri- 
butions from ao**(n — 2): 

Contribution from a ^ R (n — 2) The particle moves right and then 
left. The phase contribution stays positive. The amplitude factor 
is (1/72) 2 = 0.5. 

Contribution from a 0R L(n — 2) The particle moves right and then 
left. The phase contribution stays positive. The amplitude factor 
is (l/\/2) 2 = 0.5. 

Thus the total amplitude contribution is 0.5aoifi(n— 2)+0.5o ,rl(" — 
2) = 0.5a OL fl(n - 2) + 0.5(1 - a 0LR {n - 2)) = 0.5 

Thus we have constructive interference for both amplitudes in the transi- 
tions from set S to set . 

We need to show further that amplitudes a Li(« — 2) and clq RR {ti — 2) will not 
decrease our amplitudes for ao**(n). Let us consider clqll ij 1 — 2). Again, the 
two contributions arise from moving either RL or LR. The amplitude factor, 
as before, is 0.5. But the phase factor contributions are opposite: For RL it is 
positive, while for LR it is negative. This adds 0.5aoLL( n ~ 2) to the amplitude 
O'ORL in) and subtracts 0.5aoLL(" — 2) from the amplitude aoLfl(n). Letting 
e = 0.5a LL(n — 2), since aoLR(n— 2) and aoRL(n — 2) are two positive numbers 
summing to 1, so are aoLR(n — 2) — e and ao R L(n — 2) + e. 

A similar argument holds for the contribution from ao RR (n — 2). We have 
shown, for all even n, that a 0LR (n) and a 0RL (n) are two positive numbers 
summing to one, and hence their contribution to the probability |aoLi?(^)| 2 + 
\aoRL(n)\ 2 is at least 0.5. Note that in general, aoLi(n) and ao RR (n) will be 
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non-zero, and will also contribute to the probability at zero (though it turns out 
this contribution is small). | 

5 Conclusion 

We have defined a new kind of quantum walk with (two-step) memory, and 
investigated its properties. We see it exhibits some similarities with the clas- 
sical random walk (symmetric probability distribution, high probability at the 
origin), and other similarities with the quantum (Hadamard) walk (oscillatory 
behaviour, "tails" that propagate faster than in the classical case). We prove 
the remarkable feature of localization at the origin: in the n — > oo limit, for a 
symmetric initial state, the probability the particle is found at the origin is not 
less than 0.5. 

A referee has pointed out the work of Kendon ([§]) on decoherence in quan- 
tum walks, where probability distributions that peak at the origin are also 
obtained. However, a fundamental difference is that our peak at the origin is 
independent of the walk length, unlike the results for the decoherence case. It 
is worth examining this in more detail to see if there are other similarities in 
the results. 

Other models of quantum walks with history have been constructed (see 
[151, 14]) by using multiple coins or modifying the Hamiltonian. We find intriguing 
that the probability distribution for the 2-coin model in [H Figure 4] seems close 
in shape to our results (in e.g. Figure [2]), with again the fundamental difference 
that in our case the peak at the origin is much larger and independent of the 
walk length. Models with multiple internal states ([H|5]) have also been found 
to exhibit memory effects and localization. 
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A Appendix: Proof of Theorem 3.3 



Here we prove Theorem 3.3 We denote by Cl (respectively Cr) the number 



of clusters of Ls (respectively Rs) in a sequence of Ls and Rs representing a 
particular walk. For example, in LRLLRLLLRRL, Cl — 4 and Cr = 3. 

We examine firstly compositions of the integer Nl into Cl parts. Because of 



the phase dependence given in Lemma 3.1 we need to know how this composi- 
tion depends on iV£, the number of clusters of size one. The number of distinct 
compositions of Nl with Cl parts, and with no part of size 1 is 



Nl-Cl-1 
C L -l 



(A.l) 



(see for example [T2J page 15] ). If we want exactly one part of size 1, we take 
a composition of Nl — 1 into Cl — 1 parts, no part of size 1, and add the one 
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cluster of size 1. The number of ways we can do this is 



s> ( % L %- 1 ) ■ <"> 

In the general case, we want to add Nl clusters of size one to a composition of 
N L - Nl into C L - Nl parts, none of which is one. We can imagine having Cl 
boxes: Nl of them will be filled by clusters of size one (in Cl!/(Cj, — N^)\N^\ 
distinct ways), and the remaining Cl — Nl boxes will take a composition of 
Nl — Nl into Cl — Nl parts, without any ones. So we get 

Number of compositions / q \ / _ q _ \ \ d f / Cl 
of Nl into Cl parts with = I ATl I I _ , rl I = 

exactly Nl ones V N L J \ °£ ~ ly L - 1 J \ Nl 



(N.B. This formula does not apply in the extreme case Nl = Cl = Nl, in 
which case the number of such compositions is just 1.) For fixed values of Nl 
and Nr (so a fixed final position k), the number of walks with Cl left clusters 
(Nl of size 1) and Cr right clusters (N R of size 1) is 



(A.4) 



(Of course, Cl and Cr are not i ndep endent - they differ by at most 1.) Put ting 

and 



C L \ 


f Cr 


Nl 




N L 1 


\ n r 



in the phase factor from Lemma 3.1 the \[2 factors from Equations 2.5|2.6 



summing over Cl, Cr, Nl, N R we get the amplitude expression 

( -[\ n + N l+ N n ( Cl \ ( C R \ 

EEEE L/2 ni)[n r \. (a.5) 

Cl C R Nl Nk (Z) V N L J \ Nr J 

We now derive the specific expressions for the four possible states with final 
position k. In all cases we take as initial state | — 1, 0, 0) corresponding to a walk 
beginning LR . . . (This means in particular that Nl is at least 1, corresponding 
to the first L.) 

Final State \k + 1, k, 1) This corresponds to a walk of the form LR. . . LL. 
Since the sequence begins and ends with an L, we have Cr = Cl — 1, 
which removes the summation over Cr. In general, Nl > Cl > Nl, 
and either all three of these numbers are different or are equal. In this 
particular case, the possible values of Cl run from 2 to Nl — 1. 

Let us examine Nl and N R . For a particular value of Cl, Nl will run 
from 1 to Cl — 1 (the upper limit is not Cl because of our observation 
that the 3 numbers Nl, Cl and A?2 are either identical or all different 
from one another). N R will run from to Cr, i.e. from to Cl — !■ 



Since the expression (A.5 1 is for all walks, we need to restrict this to walks 
beginning with LR and ending with LL. For the composition of Nl Ls 
into Cl clusters, this forces the first cluster to be of size 1, and the last 
not to be of size 1. The fraction of walks whose first cluster is of size one 
is NI/Cl- Of these, the fraction that do not have a cluster of size 1 at 
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the end is (Cl — N\)/(Cl — 1)- Putting all of this together, the sum (A. 5) 
becomes 



N J2 if if (~ l)n+Nl+N * n ^ Cl ~ N ^ 

C L =2 N}=lNh=0 



(2)™/ 2 C L (C L - 1) 





C L -l 




N R 


N L J 


V N R 



(A.6) 



We separate out from the sum the limiting case: 



1. Nr = C R = N R : Using the result of Lemma 3.2 expression (A.6) 



becomes 



^ (2)«/2 N R (N R + 1) I n L l 



2N R -N L +2) 



(A.7) 



2. Nr > C R > N R : Using Lemma 3.2 the amplitude is 



N L -1 C L -1 



C L -2 



E E E 

C L =2 N L = max(l, iVjj = max(0, 
2C L - N L) 2C L -N R -2) 

Nl(C L ~Nl) ( C N \ 

Nr 



C L {C L - 1) 




(A. 



which gives us Equation 3.7 



Final State \k — l,fc,0) This corresponds to a walk of the form LR...LR. 
Since the sequence begins with an L and ends with an R, we have C R = 
Cl, which removes the summation over C R . In this particular case, the 
possible values of Cl run from 2 to Nl ■ 

Let us examine and N R . For a particular value of Cl, will run from 
1 to Cl- N r will run from 1 to C R , i.e. from 1 to Cl- We now restrict 



Expression (A.5) to walks beginning with LR and ending with LR. In 
the composition of Ls, the first cluster must be of size 1: N^/Cl is the 
fraction of walks having this property. In the corresponding composition 
of Rs, the last cluster must be of size 1: N R /C R is the fraction of walks 
having this property. Putting all of this together and applying Lemma 



3.2 the sum (A.5) becomes 



N L C L C L 

E E 

C L =2 Nl = max (l, N R = max(l, 
1C L -N L ) 2C L - N R ) 



E 



We consider the cases 
1. Nl — Cl 



(2)"/ 2 



cl 



C L \ 


/ C L 






N L J 


V N R 



(A.9) 



N£ and N R 



Cr 



Nr 



This corresponds to an alter- 
nating sequence of Ls and Rs, LRLRLR . . . LR. Clearly this path 
only exists if Nl = Nr, so the amplitude is simply 5iy L ^ R /2 n / 2 . 
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2. N L = C L = Nl and N R > C R > N R : Expression A AM becomes 




2JVr - N 



(A.IO) 



3. N L > C L > Nl and N R = C R = N R : Expression (|A.9| becomes 



«R-1 

E 

»i = max ( i ' 

2JV„ - 



(-1) 

(2)«/ 2 



Wl + WI ATI 



(A.11) 



4. N L > C L > Nl and N fl > C fl , > N R : The amplitude is 



n l -i 

E z. 

C L =2 n } = ma 



(1, ivjj = , 



(2)«/2 







Nl 




N L ) 


\N R R J 




(A.12) 



Final State \k + l,k, 0) This corresponds to a walk of the form LR...RL. 
Since the sequence begins with an L and ends with an L, we have C R — 
Cl — 1, which removes the summation over C R . In this particular case, 
the possible values of Cl run from 2 to Nl- 

For a particular value of Cl, Nl will run from 2 to Cl, while iV^ can run 
from to Cl — 1- We restrict Expression (A.5) to walks beginning with 
LR and ending with RL. This places restrictions on the composition of 
the Nl Ls, but not on the Rs. The fraction of walks that begin with a 
single L is NI/Cl- Of these, the fraction that end also in a single L is 
(Nl — 1)/(Cl — !)• Putting all of this together and applying Lemma 3.2 



the sum (A.5 1 becomes 



N L 



Cl-1 



Cl=2 N\ = max(2, nJj = max(0, 
2C L -N L ) 2C L -N s -2) 



(_l)n+yj+jyj NjjNl - 1) 
(2)"/2 C L (C L -1) 




Cl-1 
(A.13) 



We consider the cases 

1. Nl = Cl = Nl and N R = C R = N R : This corresponds to an alter- 
nating sequence of Ls and Rs, LRLRLR . . . RL. Clearly this path 
only exists if Nl — N R + 1, so the amplitude is simply 5n l -i,N R /Z n ^ 2 ■ 

2. N l = Cl = Nl and N R > C R > N R : The summations over Cl and 



Nl vanish and we get 

Nl-2 

E 



(_l)N R +N^ ( N L -l 



2Wr -J 



*(0, 
-2) 



(2)"/ 2 



No 



(A.14) 
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3. Nl > Cl > Nl and Nr = Cr = N R : The summations over Cl and 
N R vanish and we get 



(_l)iVL+JVi jy^jyi _ i) 
1 ^ (2)«/2 ^(JVh + 1) 




(A.15) 



4. N L > C L > Nl and N R > C R > N R : The amplitude becomes 



Nr 



CV-2 



E E E 

C L =2 »1 = m ax(2, N R = max(0, 
2C L - N L ) 2C L -N R -2) 



(-l)n+Nl+K N 1( N 1 _ 1) 



C L (C L - 1) 



C L \ 


C L -l 






N L 1 


V N R 



(A.16) 



Final State \k — 1, k, 1) This corresponds to a walk of the form LR . . . RR. 
Since the sequence begins with an L and ends with an R we have Cr = Cl, 
which removes the summation over Cr. Cl runs from 1 to Nl while Nl 
runs from 1 to Cl- Cr runs from 1 to Nr — 1 while N R runs from to 
Cr-1. 

For the walk to be of the required form, 

• the composition of Nl Ls must begin with a single L. This gives a 
factor of NI/C'l- 

• the composition of Nr Rs must not end with a single R. This gives 
a factor of (Cr - N R )/C R . 



Applying Lemma 3.2 the sum (A.5) becomes 



E E E* ^ n+Nl+NkN ^ CR - N ^ 



Cl = 1 = max(l, N R = max(0, 

2C L -N L ) 2C L ~N R ) 



(2)™/ 2 



cl 




We have the two following cases: 



1. N L = C L = Nl : Expression (A. 17) reduces to 



N f2 (~1) NR+Nk N L - N R { % 



N R = »ax(0, 
2N L - N R ) 



(2)"/ 2 



Nl 



Nr 
Nr 



(A.18) 



2. N L > C L > Nl : Expression (A. 171 reduces to 



(-l) n + N i+ Nl K Nl(C L - N R ) 



(2)«/2 



cl 



2C L - N L ) 2C L - N R ) 



C L \ 




Nl 


Nr 


N L J 


\ N R 



(A.19) 
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B Appendix: AXIOM code 



The following is the AXIOM source code for simulating a 40-stcp quantum walk 
with memory, and initial state | — 1,0,0). 

— In AXIOM, comments begin with two hyphens. This is a comment. 
— coin gives zero (reflection): S: |n-l,n,0> — > |n,n-l,0> 

S: |n+l,n,0> — > |n,n+l,0> 
— coin gives one (transmission): S: |n-l,n,l> — > |n,n+l,l> 

S: |n+l,n,l> — > |n,n-l,l> 
— In general, comments refer to PRECEDING line(s) of code 
P:=matrix( [[0/1 for i in 1..400] for j in 1..400]); 
— set up a 400x400 transition matrix 
for i in 1..394 I (divide (i ,4) . remainder = 0) repeat 

P(i-2,i) :=1 

P(i+5,i) :=1 

P(i-2,i+l) :=1 

P(i+5,i+l) :=-l 

P(i+4,i+2) :=1 

P(i-l,i+2) :=1 

P(i+4,i+3) :=1 

P(i-l,i+3) :=-l 

— the entries of the transition matrix: normalizatiion factor (2**0.5) left until later 
START :=matrix([ [0/1] for i in 1..400]) 
START(200,1) :=1 

— START is the initial state vector: Instead of starting at 0, we start at point 200, 

— the mid-point of the vector 

COUNT :=matrix([ [0/1] for i in 1..100]) 

— This is just for bookkeeping: COUNT will check probabilities are normalized. 



—The initial state |-1,0,0> corresponds to START(200 , 1) : =1 

—States are ordered: |-1,0,0>, |-1,0,1>, |1,0,0>, |1,0,1>, |0,1,0> etc. 



Q:=matrix( [[0.0 for i in 1..400] for j in 1..61]); 

— we will store the wavefunction amplitudes in Q 

PR0B:=matrix([[0/l for i in 1 . . 100] for j in 1..60]); 

— PR0B is the actual probabilities for each final state.... 

for j in 1..60 repeat 

— we run 60 steps of the walk 

div:=2**(j-l) 

for i in 1..400 repeat 

Q(j ,i) :=START(i,l)/(div**0.5) 

c:=0/l 

for k in 1..99 repeat 
dummy : =4*k 

PR0B(j,k):= (START(dummy, 1)**2 + START (dummy+1 , 1) **2 + START (dummy+2, 1)**2 + START (dum 
— Calculate sum of squares of the amplitudes to get probabilities. 4 amplitudes contribute 
— to the probability at each point. 

c:=c+PR0B(j ,k) 

— c should be a vector of Is, if our probabilities are normalized. 



1G 



COUNT (j , 1) :=c 
START :=P* START 
PROB 

— Finally, display the probabilities. 
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